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' Abstract. In the first part of this paper, we get new Li-Yau type gradient 

' estimates for positive solutions of heat equation on Riemmannian manifolds 

with Ricci{M) > —k, k £ M.. As applications, several parabolic Harnack 
inequalities are obtained and they lead to new estimates on heat kernels of 
' manifolds with Ricci curvature bounded from below. In the second part, we 

^Nj ' establish a Perelman type Li-Yau-Hamilton differential Harnack inequality for 

heat kernels on manifolds with Ricci{M) > —k, which generalizes a result 
I I , of L. Ni [21, 22]. As applications, we obtain new Harnack inequalities and 

C J • heat kernel estimates on general manifolds. We also obtain various entropy 

' monotonicity formulas for all compact Riemannian manifolds. 
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1. Introduction and main results 

Let (M^,g) be a complete Rimannian manifold. In the fundamental paper [18], 
Li and Yau studied the heat equation solutions 

(1) dtu = AgU 

on general Riemannian manifolds. The results in [18] has tremendous impact in 
the field of geometric analysis. One of the fundamental results is the following 
important gradient estimates for heat equations. 
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Theorem (Li-Yau [18]) Let {M,g) be a complete Riemannian manifold. Assume 
that on the ball B2R, Ric{M) > —k. Then for any a> 1, we have that 

(2) -P(^-"V)^ ^(#T + Vfci?)+^ + 



na 
2t 



Moveover, when {M, g) has nonnegative Ricci curvature, letting R-^ 00 and a ^ 1, 
(2) gives the sharp estimate (a Hamilton- J acobi inequality): 

When Ricci{M) > 0, (3) gives a clean sharp estimate. In general, on a complete 
Riemannian manifold, if Ricci{M) > —k, by letting i? ^ 00 in (3), one obtains 

|Vup Ut na^k no? 
~ - 2(a- 1) ^ 

In [7], Davies improved this estimate to 

IVuP Ut no?k no? 
f4) ■ a— < 1 . 

Let's denote the right hand side to be (^(t). Clearly, when t is big, converges to 
2 ) 

4(a-i) ■which is greater or equal to nk for any a > 1. Namely, the optimal estimate 
for t large one can get from this estimate is nfc, which can be obtained by choosing 
a = 2. For small time the dominant term of <~p{t) is By checking examples 
for heat kernels on hyperbolic spaces, one finds that when t is small, the leading 
term should be ^. This suggests one should choose a close to 1 and when time 
t is small, the sharp form is In (3) and Davies' improved estimate, if one lets 
a ^ 1, then will blow up. This phenomena suggests that there is still room to 
improve the estimate. 

It is a long time question : can one find a sharp (explicit) form for general mani- 
folds with Ricci{M) > —kl (see Problem 10.5 in book [6], page 393.) In this paper, 
we make some progress for this question. 

The first main theorem in this paper is the following local gradient estimate. 

Theorem 1.1. Let {M,g) be a complete Riemannian manifold. Let B2R be a 
geodesic hall centered at O € M. We assume Ricci{B2R) > —k with k> 0. If u is 
a positive solution of the heat equation 

{A-dt)u{x,t) = on B2rx{0,T], 

where < T < 00 and let / = In m, then we get the following Li- Yau type gradient 
estimate in Br 

(5) sup5„(|V/|2 - aft - v){x,t) < ^ + H^coth(yfc . R) + j^^^^. 



where C depends on n, a{t) = 1 + ^'"h('=^*)^™'^Mfc*) fc* ^^^^ _ |^coth(A:f) 



tanh(fet) ■ 

sinh^(fct) 

Moveover, letting i? — > 00, if Ric{M) > —k on the complete manfiold, then 

(6) |v/|2 _ (1 + sinh(feWg*)-fet )^^ ^ f[coth{kt) + l]. 
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Remark 1.1. When Ricci > 0, letting k ^ i), we recovered the celebrated sharp 
Li-Yau gradient estimates. Our estimate also explains why in the Li-Yau gradient 
estimates for general case (2), one tends to assume the blow-up parameter a > 1. 
The reason is one can view the parameter a as a function of time t, i.e. a{t) = 1 + 
sinh(fct)^cosh(fct)-fct ^ ^^^^/j indeed is greater than 1 for allt > in case of Ricci{M) > 
—k with k > 0. Since we only assume fc G R m the proof, our estimate in fact 
also works for Ricci positive case as well. However, in the positive Ricci case, 

/.\ -1 I sinh(kt) cosh(kt) — kt ^ -i 

= 1 + — sL-m — < 1- 

A linearized version of Theorem 1.1 is the following. 

Theorem 1.2. Let {M,g) be a complete Riemannian manifold. Let B2R be a 
geodesic ball centered at O G M. We assume Ricci{B2R) > —k with k>0. If u is 
a positive solution of the 

(A - dt)u{x, t) = on B2R x (0, T], 

where < T < 00 and let f = lnu, then we get the following Li- Yau type gradient 

estimate in Br 

(7) 

sup^JI V/p - af, - ^){x, t)<^ + cotHVk . R) + £2!W^!MM, 

where C is a constant depending only onn, a = 1 + ^kt and ip{t) = ^ + ^(1 + jfct). 

Moveover, letting R 00, if Ric{M) > —k on the complete manfiold, then 

(8) |v/|2 - (1 + < iL + f(l + iH). 

Remark 1.2. The global estimate (8) in Theorem 1.2 was also obtained in [1] by a 
different method. The local estimate (7) is new. Our proof seems to be simpler and 
is more of the local spirit of the classical Li-Yau's result. Moreover, the method we 

used can be extended to a matrix version. 

Remark 1.3. (8) is in the same spirit of (6) without blow-up parameter a. In 
addition, one can see from the proof that the first variation vanishes if M is an 
Einstein manifold and u satisfies the following gradient Ricci soliton equation, (a 
concept first introduced by R. Hamilton in the study of Ricci flow [12]) 

(9) ii?y-ViV,(lnu)-lffi,=0. 

In spirit, Theorem 1.1 and Theorem 1.2 are very close. The difference is the 
choice of a{t) and ^p{t). Inspecting the following series expansion of functions a{t) 
and ^p{t), one can compare these two theorems. 

^[coi\i{kt) + l\ = f^ + \nk+^{kt)-^{ktf + 0{{ktf) 

= ikt-±{ktf-Mkif+omy)- 

Indeed, functions in Theorem 1.2 are the leading terms of the expansions of func- 
tions in Theorem 1.1. Moreover, one can show by computations that 1 + '^"^^^ sinh'°(fct) *^ ~ ^* — 
1 + ffci and ^^[coth(/ct) + l] < i + |A: + ffc^t. This implies that Theorem 1.1 
yields sharper estimate than its linearized version. Theorem 1.2. 
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In [30, 31], Yau established the following gradient estimate: if Ricci{M) > —k 
with A; > 0, then 

(10) |V/p -ft< V2^^\S7f\'' + n + 2nfc + §, 

for all t > 0. Later, Bakry-Qian [1] improved the inequality to the following 

(11) iv/p -ft< v^^wTfTf + §. 

The righthand sides of (10) and (11) blow up as ^ + 0{-^) when t is small, 
while (6) and (8) give sharper estimates which has blow up order of ^ . 

In another direction, Hamilton [10] proved 



(12) \^_^2,t}H^^4,tn 

which is also sharp in the leading term for small t. But when t — > oo, the righthand 
sides of (8) and (12) will blow up, while the estimate (6) in the main theorem stays 
bounded which implies a better estimate. In regard of Li-Yau-Davies estimates 
(3)-(4) and Hamilton's estimate, one can see that the new estimate (6) works for 
both large and small time t. 

We can extend (8) and (6) to the following : under the same hypothesis of 
Theorem 1.1, the following holds 

a{t)ut + if(t)u + 2Du{y) + u\V\'^ > 0, 

for any vector field V, where a(t) and ip{t) are defined as in Theorem 1.1 and 1.2 
respectively. When fc = 0, this form of Li- Yau estimate was first pointed out in 
Hamilton's work [10]. Choosing V = 0, we get 

(13) -ait) ft < fit), for all t > 0, 

where / = Inu. One immediate application of (13) for is that (1 + ^kt)~^ e^'^^u 
is monotonic in t. When k = 0, the monotonicity of t^u is known. 



The sharp Li- Yau gradient estimate has tremendous impact in the past twenty 
years. On one hand, this gradient estimate is a differential Harnack inequality. 
Namely, it leads to a classical parabolic Harnack inequality which further yields 
powerful estimates for heat kernels on manifold with nonnegative Ricci curvature. 
There is a vast literature in studying heat kernel even before Li-Yau's work. On 
the other hand, the idea of Li- Yau leads to Hamilton's Harnack inequalities in the 
study of Ricci flow which plays a central role in his famous program. We will discuss 
more along this direction in the second part of this paper. 

Along the line of Li- Yau, we find applications of the new gradient estimates in 
deriving Harnack inequalities and new estimates on heat kernels. For example, we 
use our gradient estimates to obtain the following Harnack inequatlity. 

Theorem 1.3. IfM is a complete, noncompact Riemannian manifold with RicciiM) > 
—k. If uix,t) : M X (0,oo) K"*" is a positive solution of the heat equation on M, 
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(14) u{xi,ti) < u{x2,t2)Ai{ti,t2) ■ exp 



-{1+A2{ti,t2)) 



then for \/xi, X2 G M , < ti < t2 < oo, the following inequality holds: 

dist^{x2,x\) ^ 
4(i2 - ii) 

where xi,X2 S M, < ti < t2 < oo, dist{xi,X2) is the distance between x\ and X2, 

A _ f £^*2-2kt2-l\^ „„J A (+ + \ — t2 COth(fct2)-tlCOth(fcti) 

M - [ e^kH_2ku-i ) , and A2{ti,t2) • 



Remetrk 1 



,4. Easy to see lim Ai{ti,t2) = f — ) ^ , lim A2{ti,t2) = 0. 



Similarly, the linearized gradient estimate also yields a corresponding Harnack 
inequality. 

Theorem 1.4. If M is a complete, noncompact Riemannian manifold with Ricci{M) > 
—k. If u{x, t) : M X (0, oo) IR+ is a positive solution of the heat equation on M, 
then for Vxi , X2 S M, < ti < t2 < oo, the following inequality holds: 
(15) 



^{x^,t,)<u{X2,t2){^)'{]^±;^) ^exp( 



/ dist (x2,xi) 1,, ,^ n,. 



where x\,X2 € M, Q <t\ <t2 <oo, dist{xi,X2) is the distance between x\ and X2- 



One should compare the above theorems with others when t is large or small and 
when d{xi,X2) is large or small. When fc = 0, they reduce to the classical result. 

As standard, we find a lower bound of the heat kernel as well by using the 
Harnack inequality. 

Theorem 1.5. Let M be a complete (or compact with convex boundary) Riemann- 
ian manifold possibly with Ricci{M) > —k. Let H{x,y,t) be the (Neumann) heat 
kernel. Then 
(16) 



H{x,y,t)> (47rt)"t2-T 
and 



(2kt) ^ 



(e2 



-2/st-l) 4 



• exp 



kt coth{kt) — 
ki 



H{x, y, t) > (4^<)-t exp [ - %)^(1 + ^kt) - ^kt\ , 

for all x,y G M and t > 0. 

One should compare this theorem with CoroUay 2.3 in [10]. 

Remark 1.5. By going through Li-Yau's paper carefully, one can get similar re- 
sults on the estimates of Greene's function, lower bounds of Dirichelet or Neumann 
eigenvalues, betti numbers, etc. The new contribution will be that explicit depen- 
dence of various constants can be established. 



In the second part of this paper, wc will discuss the relation between Li-Yau 
type gradient estimate, Hamilton's gradient estimate, and Perelman type differen- 
tial Harnack inequality. Motivated by Li-Yau's fundamental work, Hamilton proved 
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the following gradient estimate. 

Theorem (R. Hamilton [10]) Let {M,g) be a closed Riemannian manifold with 
Ricci{M) > —k. Let u{x,t) be the positive solution to the heat equation. Assume 
u < A, then 

(17) t\Wu\'^ <{l + 2kt)u^ln{^). 

Ground breaking progress in the study of Ricci flow and complete proof of 
Poincare conjecture was made by G. Perelman in 2002-2003. Some important tools 
which enable Perelman to make the breakthrough were related to Li-Yau [18], and 
Hamilton's earlier work [10, 12, 11], (see also [4]). More specifically, Hamilton 
systematically studied the differential Harnack inequalities in Ricci flow along the 
line of Li-Yau. Perelman discovered a new sharp differential Harnack inequality for 
Ricci flow which plays a crucial role in his work. One new feature of Perelman's 
work is that no curvature assumption is assumed. Moreover, Perelman's differential 
Harnack is modelled on shrinking Ricci soliton and works for all dimension. 

A natural question is whether Perelman's new discovery could shed some lights 
on the results for linear heat equations. Indeed, one could find highly similarities 
between backward conjugate heat equation along Ricci flow and heat equation so- 
lutions on static Riemmannian manifolds. In [21], one of the main results is the 
following analogue of Perelman's differential Harnack inequality for heat kernels. 

Theorem (L. Ni [21, 22]) Let {M,g) be a closed Riemannian manifold with 
nonnegative Ricci curvature. Let u{x, t) = H{x, t; y, o) be the positive heat kernel. 

Then 

(18) t(2A/-lV/l2) + /-n<0, 
— / 

where u = — — w . 

(47rt) 2 

Easy to sec, this type of differential Harnack quantity t(2A/ — | V/|^) + / — n is 
a hybrid of Li-Yau's estimate on [V/]^ — Q;(A/ -t- [V/]^) and Hamilton's estimate 
on I V/]^ + (t + 2A;)/. As we have seen in section 3, from Li-Yau type gradient esti- 
mate, one could get a Hamilton- Jacobi inequality which leads to the generalization 
to a classical parabolic Harnack inequality of Moser. This powerful method was 
started by Li and Yau. Hamilton extended this method further for heat equations. 
Moreover, he established similar estimates in the study of Ricci flow. This method 
now is generally referred as Li-Yau-Hamilton estimate (LYH) (cf. [23]). 

In regard of the nice curvature free feature of Perelman's LYH type differential 
Harnack inequality under Ricci flow, and our new discovery of Li-Yau gradient 
estimate, one may ask : can one find a Perelman type of differential Harnack 
inequality for heat kernels on any closed Riemannian manifolds? We answer this 
question affirmatively. The following is the second main theorem in our paper. 

Theorem 1.6. Suppose M" is a closed manifold. Let u be the positive heat kernel 
and k>0 is any constant satisfying Rij{x) > —kgij for all x G M", then 

(19) V := [tAf + t{l + kt){Af - ]V/]2) + / - n(l + ^ktf]u < 0, 
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for all t > with u = — — n-. When k = 0, this theorem is due to L. Ni. 
Moreover, 



(20) (^ - A)v = -2t| V,V,/ - + '^)9^J - 2t(i?y + kg^j)f,fjU. 

Remeirk 1.6. The right handside of the evolution equation of v vanishes if the 
manifold is Einstein and f satisfies a gradient Ricci soliton equation. See discus- 
sions in section 6. 



Remark 1.7. The proof of the above Perelman type LYH differential Harnack in- 
equality is independent of the Li- Yau type estimate. In fact, besides the discovery of 
this differential Harnack quantity for general manifolds, we have overcome a tech- 
nical difficulty in the proof, which in the Ricci{M) > case, see [21, 22], previously 
one has to appeal to Perelman' s reduced distance function, see Remark 4-1- 



Point-wise differential Harnack inequalities and monotonicity formulas for en- 
tropy functionals arc closely related. Usually, a point-wise differential Harnack 
quantity easily yields a monotonicity formula for the related functional. But re- 
versely, it is more difficult to find the corresponding differential Harnack quantity 
from a functional monotonicity. In this paper, we will analyze this relation and 
give various different new entropy monotonicity formulas for heat solution. In fact, 
this served as one of the motivations of this paper. 



We introduce the following Li- Yau type entropy formula Wly and Perelman 
type entropy formula Wp, which were formulated from our new point- wise Li- 
Yau differential Harnack quantity and the new Perelman type differential Harnack 
quantity respectively. 



y\^LY{u,t) 



or 



(21) 



/ 

I 

JM 



smh^ {kt) 



^[coth(fct) 



1 



udfj, 



\ , n nk 1 , , 

Alnu + - + —{l + -kt) udn, 



1, and 



^dfi, 



u in Wp. 



Wp(/,r)= / (r|V/|2 + /-n(l + ifcT)^ 

where u is a positive solution of the heat equation, ^ = 

We shall prove that these entropy functionals are nonpositive and monotonically 
nonincreasing for all t > on manifolds with Ricci{M) > —k. Moreover, the first 
variation vanishes if and only if the manifold is Einstein and In u satisfies a gradient 
Ricci soliton equation. 



We also discuss estimates for Nash entropy on a closed manifold with Ricci{M) > 
-k. 



Remark 1.8. As in Li-Yau, one could extend all the results in this paper to heat 
equations with potentials. In particular, one can obtain better Harnack inequality. 
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and lower, upper bounds for the fundamental solution. This was treated by the 
authors in a separated paper [15]. 

This paper is organized as following. In section 2, we prove the generalized Li- 
Yau gradient estimates for manifolds with Ricci curvature bounded from below. 
In section 3, we discuss applications of the gradient estimates and obtain classi- 
cal Harnack inequalities and estimates for heat kernels. In section 4, we prove 
the Perelman type Li-Yau-Hamilton Theorem 1.6. In section 5, as applications of 
Theorem 1.6, we derive another parabolic Harnack inequality for heat kernels. In 
section 6, we discuss various entropy formulas with monotonicity for heat equations. 

In this paper, we will use Einstein convention. 

2. Li-Yau type gradient estimates on general manifolds 

We start with the following lemmas. 

Lemma 2.1. Let M" he a Riemannian manifold, and u{x,t) is a positive solution 
of the heat equation. Let a{t) and (f{t) be functions depending on t and F = 
|V/p-a/t-^. Iff = \nu, then 

ft = A/+|V/|2 
iA-dt)ft = -2V/V/* 
^^^^ (A-5*)|V/|2 = 2|/,,f -2V|V/|2V/ + 2i^,,/,/,■ 
(A-a^)F = 2\fij\^-2VFVf + 2Rijfifj + a'ft + ^' 

Proof. Prom Ut = Au and f = \nu, we get ft = A/ + |V/p. Hence, 

(A-at)(A/ + |V/n = A(A/ + |V/|2) - A/i - 2V/tV/ 

= -2V/tV/. 

The second identity follows from Bochner formula. 

(A-5t)|V/|2 = A|V/|2-2V,AV/ 

= 2|/i,f + 2VA/V/ + 2R,,fJ, - 2V/tV/ 
= 2|/yf - 2V|V/|2V/ + 2Rijfifj. 

The last identity in (22) follows from the definition of F and the first three identities. 

□ 

Lemma 2.2. Let M" be a Riemannian manifold. Suppose u, f, and F are defined 
the same as in Lemma 2.1. By choosing different sets of a{t) and (p{t), we have 
the following. 

(1) Ifa{t) = 1 + ffci and ip{t) = § + f A; + ^kH, then 

(23) (A - dt)F = 2\fij + l^gij + ^gi^ |2 - 2VFV/ + f F + 2{Rij + kgij)fifj 

(2) // a{t) = 1 + '^"^^^tinfffc?^"'* ^(*) = [coth(A;t) + l] , then 

(24) 

(A - dt)F = 2\f,j + f ,9yf - 2VFV/ + 2fc coth(fcf)i;^ + 2{R,, + kg,,)fjj 

Remark 2.1. By looking at the series expansion, one can see that the first set of 
a{t) and tp{t) are the linear approximation of the second set. Moreover, one can 
show that 1 + '"'''^tiSfc?^"''* ^ 1 + i't* f[coth{kt) + 1] < § + f fc + ^kH. 
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^<p-2k= 2kcoth{kt) 

— <p—a' 

a = 



\(p-2k 



Proof. We only prove (24). The proof of (23) is similar. From lemma 2.1, we get 
(25) 

(A - dt)F = 2|/y |2 - 2VFVf + 2Rijfifj + a' ft + ^' 
= 2|/,,- + ^3i,f -^¥'A/-^-2VFV/ 

-2fc|V/|2 + 2{Ri, + kgij)fifj + a' ft + if' 
= 2\fij + !^gij\^ - 2VFV/ + 2{Rij + kgij)fifj 

-2fc|V/p ^ ^Mft - |V/P) - ^ + a' ft + ^' 
= 2\fij + ^gij\^ - 2VFV/ + 2{Rij + kgij)fifj 

+(^^ - 2fc)(|V/p - 1^/, - ^) - ^ + ^' + (1^ - 2k)^ 

Direct computations yields that, a{t) = i+ '=°^MfcJ)^BmMfc*)-fct ^nd (pit) = [coth(fci)- 
l] , satisfy the following system 

(26) 

Plug into (25), we get 

(A - dt)F = \fij + ^gij]'' - 2VFV/ + 2fccoth(fct)F + 2(i?,,- + kgij)fifj, 
which completes the proof. □ 
We also have 

Proposition 2.1. Let M" be a Riemannian manifold. Suppose u is a positive heat 
solution, / = — Inw, and F = |V/p + a{t)fi — (p{t). By choosing different sets of 
a(t) and (p(t), we have the following indentities respectively. 

(1) Ifa{t) = l + lkt and fit) = § + ffc + %kH, then 

(27) (A - dt)eFu = 2t^\f^, - - |5«,f w + 2t^iRij + kg,j)fj, 

(2) // a{t) = 1 + ^{t) = f [coth(fct) + 1] , then 

(28) 

(A - dt) smh'{kt)Fu = 2 smh''{kt)\fij + ^gij + 2 smh\kt){Rij + kgij)fifj. 

Proof. Wc will only prove the first identity. The second one is similar. Apply 
Lemma 2.2 to m = e^-' , wc get 

(A - dt)F = 2\- fij + l^i,- + ^5i,f + 2VFV/ +jF + 2{Rij + kgij)fifj 

Using Ut = Am and direct computations, we can show that 

(A - dt)t''Fu = 2t''\fij - ^ft, - '^gijfu + 2t\Rij + kgij)fifj. 

□ 
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Consequently, the following estimates on closed Riemannian manifolds hold. 

Theorem 2.1. Let M" he a compact Riemannian manifold possibly with boundary 
and with Ricci{M) > —k. Suppose u{x,t) is a positive solution of the heat equation. 
If dM ^ 0, assume that DM is convex, and u{x,t) satisfies the Neumann boundary 
condition 

du 

— =0 on dMx(0,oo), 
ou 

where |^ denotes the outer normal of dM. If we let f = \nu, then we have 

(29) |v/|2 - (1 + SiHMMlE;^*)^)^^ < i^[coth(fct) + l]. 
On the other hand, the following linearized version is also true. 

(30) |V/|2 - (1 + |fct)/i < IL + |fc + |fc2^. 

Proof. The proof is by the standard parabolic Maximum Principle (cf. [26]). We 
will skip the details. □ 

To prove the main theorem Theorem 1.1, we need the following technical lemma. 

Lemma 2.3. For x > 0, the following is true, 

k 

4(a - 1)g9 - na^ (2fccoth(a;) r^^-^-r^) < 0, 

^ cosh(a;) sinh(x) ' 

where a{x) = 1 + ""^^^^ii^^''^'" and ^{x) = f ( coth(x) + l) . 
Proof. Equivalently, we need to prove 

2(a - l)(Sgg^ + I) - ^"uK^^u\ < 

2(a- 1) • (e^^ + 1) -a2(e2^ + e-2^) < 

2(a-l)-(p>' + e2^)-a2(e4^ + l) < 
-e^^(a-l)2-e'^^ + 2(a-l)e2^-a2 < 

-e4^(a - 1)2 - [e2^ _ (a - 1)] % (a - 1)^ - < 

Since a = 1 + '"'"^^J.^J^T^'j)"'^"'" , we get < a - 1 < a. Hence the last inequality is 
true, which finishes the proof of this lemma. □ 

We now prove the main theorem, Theorem 1.1. 



Proof (Proof of Theorem 1.1) Let's denote F = |V/|2 - a{t)ft - (p{t) and 
G = (3{t)F, whore a{t) and ip{t) are defined as in the main theorem and (3{t) is a 
positive function of t to be determined. Applying (1) in Lemma 2.2 to G = l3{t)F, 
we obtain the following 

(31) 

(A - dt)G = 2(3\fij + ^gij\^ - 2VGV/ + G{2kcoth{kt) - ^) + 2/3(i?y - k)fifj 
> 2(3\fij + ^gij\^ - 2VGV/ + G{2kcoth{kt) - f ). 
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Recall that one can construct a cut-off function (j) as in the proof of Theorem 4.2 
in [26], which satisfies supp (j) C B2R and = 1. Moreover, 

(32) 



A(?i > -%{l + ^/k-Rcot\i{^/k-R)^, 



where C depends only on n. We want to apply maximum principle to (j)G on B2R x 
[0,T]. If (/)G attains its maximum at (a;o,io) <= -B2_r, x [0,r], then ((i)G)(.To, to) > 
without loss of generality. So xq € iJ2_R5 io > 0, and by maximum principle, at 

= V{(j}G) = GV(j) + (ING 

(33) A(^G) < 

at(<^G) = <^Gt>0. 

Notice that for an n x n matrix A, we have \A\'^ > ^ (trA^ , and 

(34) tr(/i,- + ^5,,) = A/ + = -1 [| + (a - 1)(|V/|2 - <p)" 

In the sequel, all computations will be at the maximal point (xq, to) and we will 
frequently use property (33) whenever necessary. Applying (31), we have 

0> (A-at)(</.G) = GA</.-2GM + ^(A-a,)G 

> g(a</, - 2^) - 2GV<^V/ + 2cj>l3\fi^ + ^gi.f + <^G • {2k coth(/ci) - f ) 
Multiplying by (j), and applying (34), we have 

0> (0G)(A(/)-2l^ + 2fccoth(H)(^- ^) - 2((?!.G)|V0||V/| 



.20^/3 
noc^ 



f + (a-l)(|V/p-(^) ' 



= (,^G)(A,^-2M + 2fccoth(fct)</.-f) - 2(0G)|V</.||V/| 

+ ^,G^ + ^^^^^(|V/P - ^)^ + m . ^^(|V/p - ^) 
= (0G)(a</, - 2^ + 2fccoth(H)<^ - f - + 

+ ^^^^^(|V/p - + m . (^^|V/|2 - 2|V.^||V/|) 

Applying inequality ax"^ — bx > —j^, {a > 0), to the last term, and also drop 
the second last term which is nonnegative, we get 

> (^G)(a0 - 2^ + 2fccoth(fci)0 - f - - ^g^) + ^,{<PGr. 
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(35) 



Since (0G)(xo,io) > 0, we get 

(0G)(a;o,<o) < f 4(a - 1)^ - W (2fc coth(H) - 



^f^^4, + 4<f.(Q-l) 



where the righthand side is evaluated at {xo,to) which depends on the function. 
Hence on Br x [0,T], applying estimates (32) on (f), we have 

G{x,t)< {<pG){xo,to) 



(36) < I 



4(a - l)(p - na^ (2fc coth(fct) - ^ 



/3 ^ 



+ ^ (f + ^ COth(V^ • i?) + 4^^) . 



Next, wc choose I3{t) = tanh(fci). Hence, ^ = sinh(fct) rosh(fct) • Denote x = kt 
and applying Lemma 2.3, wc have 

(37) 4(a - 1)(/? - na^ {2k coth{kt) - ^) <(). 

On the other hand, by definitions, P(t) = tanh(fct), a{t) — 1 = ^ 
as t ^ 0, we have 2: as t ^ oo, we have 1. This implies 

' a—l ' ' a—1 ^ 

(38) ^ < C, 

a—l 

where C is a universal constant. 

Recall that all the computations arc at {xo,to) and (xo,to) is the maximum 
point, to <T and P{t) = tanh(A:f) is non-decreasing. Plug (37) and (38) into (36), 
we get 

i^G){x,T) < {(t>G){xo,to) < n2!(*|M(*o) + c^coth(x/fc . i?)) + ^ 
< /J(T)(^ + 2^coth(^/fc•i^)) +2^, 

where in the last inequality, we have used the fact that a{t) is uniformly bounded 
over (0, oo). But ^=1 on Br, hence, from G = pF, we have 

supj,^F{x,T)< ^ + n^coth(^/fc•fi)+ ^,,^:^(,^) , 

since T is arbitrary, the theorem has been proved. 

□ 

Similarly, choosing a different set of a{t) and <f{t), one can prove the linearized 
local version. Theorem 1.2. 

Proof. (Proof of Theorem 1.2) The proof is similar to the proof of Theorem 1.1. 

The diff'erence is the choices of if{t) and a{t). Denote F = |V/p — a{t)ft — ^{t) 
and G = (3F, where a{t) = 1 + |A;i and ip{t) = H + ^(1 + \kt) are defined as in 
Theorem 1.2. Applying (1) of Lemma 2.2 to G = fi^, we obtain the following 
(39) 

(A - dt)G = 2p\fij + (i + |).9„ |2 - 2VGV/ + f f - ^) G + 2f3{Rij - k)fifj 
> 2/3|/,, + (i + |)5i,p-2VGV/ + 
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Construct the cut-off function (j) as before, which satisfies supp (p C B2R and 
(f)\BB, = 1. Moreover, 



\W4 



< 



fl2 



(40) 



> -§,(l + \/k-Rcoth{\/k-R)y 



Apply maximum principle to (f>G on B2R x [0,T]. If attains its maximum 

at {xo,to) G x [0,T], then {(f)G){xo,ta) > without loss of generality. So 
xq € B2R1 to > 0, and by maximum principle, at (xq, to) 



(41) 



= V(^G) = GV^ + 0VG 
A(0G) < 
dti<l)G) = <j)Gt>0. 



Using the matrix inequality > -^(trA^ , we have 



(42) 



A/+f(i+.) = -i 
2 t a 



^ + ia-l)\Vf\'-{^ + lnkH) 



In the sequel, all computations will be at the maximal point {xo,to)- Applying 

(39), wc have 

0> {A-dt)m = GAct>-2Gi^ + <P{A-dt)G 



> G 



,|V0|^ 



. 2GVW + '^Mfv + iTt + l)9vf + (f - 'j)^G 



Multiplying by 0, and applying (42), wc have 
> (0G)(a<^- 2M + (I _ |)<^) _ 2(,^G)|V<^||V/| 



+ 



f + (a-l)|V/p-(f + infc2t) 
= (0G) (a,/, - 2^ + (f - !)</>) - 2(</,G)|V</.||V/| 



G2 



20-' /3 

■ 2~ 



(a-l)|V/|2-(f + Wi) 



(0G) 



4^ 



(a-l)|V/|2-(f + inft2f) 



= (<^G)(a0 - 2^ + (f - f )</, - ^(f + lnkH)cp) + ^G^ 



+ 



2^/3 



(a-l)|V/p-(f + infc2t) +(,^G).( 



40(a-l) 



|V/p-2|V<^||V/|) 



Applying inequality ax"^ — bx > —j^, {a > 0), to the last term, and also drop 
the second last term which is nonnegative, we get 



> 



m[A4>- 



+ (f-f 



2 /2/s , 1 p.w _ na^\V4>\^ \ , 2 (J.n\2 

Ks- + 3'^ 40(a-i) ; + ;iz^5^v9'^^ • 
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Since {(pG){xo,to) > 0, we get 

/.^N/ X na^/jr /3' 2. 2,2k 1,2x1. na^f3 / „|V(/)P na^ly^pN, 

(^G)(.o, to) < [(^ - i ) + ^(y + 3^^*)] ^ + ( - + 2^ + ^^(^) 

where the righthand side is evaluated at {xo, to) which depends on the function. 

Hence on Br x [0,T], applying estimates (32) on </>, we have 
(43) 

G{x,t)< {(l)G){xo,to) 



< nl 

- 2t 



ta^f - f ) + 2C-M + 1 fc2,2)l ^+n^f^ + c^ eoth(v^ .R) + ^^ 



Next, we choose f3(t) — tanh(fcf). Then, — sinh(fct) cosh(fcf ) ' Denote x — kt and 
recall a{t) = 1 + It is not hard to show that for a; > 0, (see the comments 
after the proof of this theorem) , 

(44) ( sinh(.rcosh(x) - 2)(1 + 1^ + 1^') + (1^ + 1^')) < 0, 

this yields that 

On the other hand, by definitions, as f — > 0, P{t) = 0{t) and a{t) — 1 = 
This implies 

(45) ^ c-, 

a — I 

where C is a constant. 

Recall that all the computations are at {xo,to) and {xo,to) is the maximum 
point, to < T, a{t) = 1 + ^kt and P{t) = tanh(H) are non- decreasing. Plug (44) 

and (45) into (43), we get 



{^G){x, T) < (0G)(xo, to) < {w + ^ coth(Vfc • R)) + ^li!^ 

< "«'^(^/(^) + coth(Vfc • R)) + 

Since ^ = 1 on Br, we obtain the following gradient estimate in Br 

snps.F{x,T)< ^ + ^^coth(^/fc.i^) + ^^^. 

Since T is arbitrary, we proved the first part of the theorem. 

If the manifold is complete, for any fixed T > 0, letting i? — > oo, we get 

F{x,T)< 0. 

Since T is arbitrary, equivalently, we obtain the global estimate (8). □ 

One way of proving (44) is to prove an equivalent inequality as follows, 

(46) I{x) := (e2-- - e-2-)(l + fa; + ^x^) - 2x{l + fa; + > 0, 

where function I[x) is a real analytic function and all the coefficients of its Taylor 
expansion are positive. 
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3. Harnack inequality and heat kernel estimates 

Along the line of Li-Yau, as an application of the gradient estimates in section 
2, we can establish Harnack inequalities for positive solutions of the heat equation 
and deduce lower and upper bounds for the heat kernel. 

We prove the Harnack inequality for noncompact manifold first. 

Proof. (Proof of Theorem 1.3 and Theorem 1.4) Let / — \nu, then com- 
bining estimate (6) and the heat equation, we have the following Hamilton-Jacobi 
inequality, 

/*>-^Mt)-|v/p). 

a 

Let 7 be a shortest geodesic joining xi and X2, 7 : [0,1] M, 7(0) = X2 
7(1) = xi. Define a curve -q in M x (0, 00), tj : [0,1] ^ M x (0, 00) by 77(5) = 
(7(s),(l - s)t2 + sti). We have r]{0) = {x2,t2), = {xi,ti). If p = d{xi,X2), 
then I7I = p. We have 

f{x,,h)-f{x2,t2) = [' ^f{v{s))ds 



< 



ds" 

[{^,Vf)-{t2-h)ft)ds 

I' (p|v/l + ^^(¥'W-|v/p))ds, 



where t = (1 — s)t2 + sti. 

The integrand is a quadratic polynomial in |V/|, whose maximum value is 



ap"^ t2 - t 



A{t2-h) a 



, "-a ~ "-i /,\ 



Therefore, we obtain 



f{xiM)-f{x2,t2) < f^(j^^^^a{t) + ^-^^{t))ds 



4{t2-t 



dt 



_ jj_ I fctcoth(fct)-l \ I 7t 1 smh(2fct) -2fct+cosh(2fct) - 1 

~ [4(t2-ti)M'' + fc j-r^"^ 2k 

where in the second identity we have used t{s = 0) = ^2, t{s — 1) — ti, dt — — {t2 — 
ti)ds, and we have chosen a(t) = 1 + '^'"^''^g*^„^T(fe^)'''~*'^ and if (t) = ^ [ coth(A;t) + 1] . 
Taking the exponential of both sides and flip the quotient, we finish the proof of 
Theorem 1.3. Similarly, choosing a = 1 + and ip{t) = ^ + ^(1 + ^kt), we 
prove Theorem 1.4. 

□ 
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Theorem 1.3 extends Li-Yau Harnack estimate in several ways. When fc = 0, we 
recover Li-Yau's theorem of the case Ricci{M) > 0. When Ricci{M) > ~k and 
A; > 0, this theorem finds the expHcit form for Li-Yau's original theorem without 
parameters and it improves previous estimates. When Ricci{M) > 0, this formula 
is new. In the last case, the Harnack is true only for short time. 

For a positive solution on compact manifolds, the Harnack inequality given by 
Theorem 1.3 also holds. 



Theorem 3.1. IfM is a compact Riem.annian manifold possibly with boundary 
and with Ricci{M) > —k. If dM ^ 0, we assume dM is convex. u{x,t) is a 
positive solution of the heat equation on M, and |f = if dM ^ 0. 



(47) u{xi,ti) < u{x2,t2)Ai{ti,t2) ■ exp 



dist^ (x2,xi) 



{1+A2{ti,t2)) 



4(i2 - ii) 

where xi,X2 & M, < ti < t2 < oo, dist{xi,X2) is the distance between xi and X2, 

A _ ( £^f£_;2fet2^\ T „^ . A /, , \ _ t2 COth(fct2)-tlCOth(fctl) 

Ai - [ ^■.kt,_2kt,-i ) , and A2{ti,t2} • 



Proof. The proof is similar to the proof of Theorem 1.3. We skip the details. 



□ 



Applying the linearized gradient estimate Theorem 1.2, we obtain another Har- 
nack inequality on compact manifolds as follows. The proof is similar which we will 
also skip. 

Theorem 3.2. If M is a compact Riemannian manifold possibly with boundary 
and with Ricci{M) > —k. If dM ^ 0, we assume dM is convex. u{x,t) is a 
positive solution of the heat equation on M, and = if dM ^ 0. Then 



(48) u{xi,ti) < ■u(x2,t2)(^) ' -exp 



dist^(x2,xA , 1,, ,x n,, 
— , ^ ' (l + -fc t2+^i ) + fe-^i 
4(^2 -ti) ^ 3 ^ " 4 ^ ^' 



where xi,X2 S M, < ti < t2 < oo, dist{xi,X2) is the distance between x\ and X2. 

It is well-known that Harnack inequality leads to lower bounds on the heat ker- 
nel. Applying the Harnack estimates, we get 

Theorem 1.5. Let M be a complete (or compact with convex boundary) Riemann- 
ian manifold possibly with Ricci{M) > —k. Let II{x,y,t) be the (Neumann) heat 
kernel. Then 
(49) 



4t 



Hix,y,t)> (4.t)-t2-f ^^J^-^^, .exp 
and 

H{x, y, t) > (47rt)- 1 exp [ - + ^kt) - ^kt 



llf-t j_ ktcoth(kt}^\ 



for all x,y € M and t > 0. 
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Proof. Since the proofs are similar, we will only prove the first inequality for com- 
plete noncompact manifolds. Applying the Harnack inequality (14) to the function 

u{y,t) = H{x,y,t), 

we obtain 



H{x,x,s)< H{x,y,t + s)[ ^''''^;L-_Ts^ 



■ exp 



d\x2 



it 



£1) ^ 



1 



(t+s) coth(fc(t+s))-s coth(fcs) 



for all s > and t > 0. By local calculations one gets 

limiA-Ks)^ H(x,x,s) = 1, 

s— >0 

SO, multiplying by lims_>o(47rs)t , we get 

(47rs)* -fr(x, X, s) < {A7Tt)^H{x, y,t + s) 



s 2 



(e2fe(«+«) _2A;(t + s)-l)f 



■ exp 



4t 



(t+s) coth{k{t+s)) — s coth(fes 
i 



t2 



Let s — > 0, we obtain 

. (e2'=*-2fci-l)^ 



1 < {A'Kt)^H{x,y,t)2-' 



This completes the proof. 



{2kty- 



■ exp 



(P{x2tXi) ( ktcoth{kt) — 1' 



4i 



1 + 



kt 



□ 



As a direct corollary of Theorem 1.5, we obtain the following estimate for 
H{x, X, t). 

Corollary 3.1. Under the same assumption as in Theorem 1.5, we have 

H(x,x,t)> (47rt)-t2-t ^ „ 

and 



H{x,x,t) > (47rt) 2 exp 



}kt 



Remark 3.1. Results in this section hold for compact Riemannian manifold with 
or without boundary condition. If the boundary is nonempty, we need to assume 
the boundary is convex. 

4. A Perelman type differential Harnack inequality 

We will devote this section to the proof of Theorem 1.6. Theorem 1.6 may be 
the closet one to the differential Harnack inequality discovered by Perelman along 
Ricci flow. It's worthwhile to note that in this section we will follow a different 
notation which was used before by Yau in [31] and Perelman in [24]. Namely, we 



{A-Kt) 2 



IT for positive heat 



always assume u = e for positive heat solutions and u ■ 

kernels. We start with a lemma. Let 

X{f,t)= tAf + f + ip{t), with (^(t) = -f ln(47ri)-n(l + ifct)2 
y{f,t)= /t = A/-|V/|2, 

z{f,t)= x + t{i + kt)y. 
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Easy to see, j'^^yudv = Jj^ ftudji = jj^^^f — \V f\'^)udv = 0. Hence, a multiple 
of Jj^ yudv with a time function will not affect the entropy functional. 

Lemma 4.1. Let M he a complete Riemannian manifold. Let u be a positive 
solution to the heat equation with u = . Then 

(1) (f -A)A' = -2t|ViV,/-(^ + |)5,,f -2VA'V/-(2fct+l)(A/-|V/|2)- 

2t{Rij + kgij)fifj 

(2) - A)3; = -2vywf 

(3) - A)Z = -2t\ViVjf - + |)5i,f - 2V2Vf - 2t{Rij + k9ij)fjj 

Proof. By direct computations with the help of Bochner formula, (c.f. Lemma 2.2 
in [21]), one can prove equality (1), 

(i-A)A' 

= -2iA|V/|2 + (A/-|V/|2) + (^' 

= -2t\fij\^ - 2tVAfVf - 2tRijfifj + (A/ - \Vf\^) + <p' 

= -2i|/i,f - 2VXVf + 2|V/p + 2fci|V/|2 + (A/ - |V/n - 2t{Rij + kgij)fifj + ^' 
= -2t|/y - + |),g,,f - 2VXWf - (2kt + 1)(A/ - jV/p) - 2t{R,, + kg,,)fifj. 

Equality (2) is again by Bochner formula and it has already been proved in section 
2 Lemma 2.1. Inequality (3) is immediate from equality (1) and equality (2), 

- A)2 = (f - A)X + t(l + - A)3; + (1 + 2kt)y 

= -2t\\7iVjf - ( + ^)gij |2 - 2V^V/ - 2t{Rij + kgij)fifj 

□ 

Now, we recall Theorem 1.6. 

Theorem 1.6 Suppose M" is a closed manifold. Let u he the positive heat kernel 
and k>0 is any constant satisfying Rij{x) > —kgij for all x G M", then 

(50) V := [tAf + t{l + kt){Af - |V/n + / - n(l + \ktf]u < 0, 

— / 

for allt>0 with u = — — w . 

•' (47rt)7 

We first derive the evolution equation equation of v. 
Proposition 4.1. Let u, f, andv he defined as in Theorem 1.6. Then 

^^^^ ^§t ~ ^^"^ ^ -2i| ViV,-/ - + ^)gij \\ - 2t{Rij + kgij)fifjU. 

Proof. Recall in Theorem 1.6, u = -^j-^. One can use the change of variable idea, 
or simply observe that / + f ln(47rt) satisfies the assumption of Lemma 4.1. Easy 
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to see v{f, t) — Z(f + ^ ln(47rt), t)u. Hence, by Lemma 4.1 



A){Zu) 



(52) 



= {^-A)Zu-2VZVu + Z{§-A)u 
= A)Zu + 2uVZWf 

= -2t\ViVjf - + - 2t(i?y + kgij)fifju, 



where we have used the fact that (/ + f ln(47rt))j = fi 



□ 



The way we find the point- wise differential Harnack quantity can be used to find 
similar quantities for other geometric evolution equations, e.g., along Ricci flow 
equation. See an application in [16]. 

One can rewrite t; as t; = [i(2A/ - |V/p) + f - n]u + kt^Au - nkt{l - jt)u. 
Clearly, the first term is the one discussed by Ni for the Ric > case and the second 
term is divergence free. The last term contributes to the Ricci curvature term from 
Bochner formula. On the other hand, the difference between v and the integrand 
of Wp, cf. (60), is a divergence term {t + kt^){Af - |V/ p)u, = (t + kt^)Au. This 
term is crucial for finding the pointwise differential inequality. The evolution of v 
also yields a proof to the monotonicity of Wp in section 6. 

When A: = 0, our theorem reduces to the result of L. Ni in [21] and [22]. The 
reason we discuss the case of fc > is because we will make essential use of the heat 
kernel comparison with hyperbolic space and Euclidean space, sec Theorem A. 3. 
in the Appendix. Going through more complicated computations, one could deal 
with Ricci{M) > case as well. Using ideas from this section, we also established 
a direct proof for Perelman's differential Harnack inequality along Ricci flow in [16]. 



Remeirk 4.1. Our proof of this differential inequality is different from the one of 
L. Ni for Ric > case in [21] and [22]. The main simplifi,caMon is we do not need 
to use the gradient estimates for all positive solutions to the heat equation, and 
also other techniques such as the nontrivial reduced distance function introduced by 
Perelman in the study of Ricci flow. 

Proof (Proof of Theorem 1.6:) For any to > 0, let h be any positive function. 
We solve the backward heat equation {dt — A)h{y, to—t)=0 starting from to with 
initial data h. We then have that. 
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dt I vhdfj. 

JM 



{hfV + hvt)dn 



M 



(53) 



/ (-A/iw + hl^v)dn -2t f I ViVj/ -{^^ + ^)gij\^hudiJL 

JM JM 2 

-2t / {Rij + kgij)fifjuhdn 

JM 

f \ I k ,2 f 

-2t / I ViVj/ - {-^ + -^)9ij I hudu -2t {Rij + kgij)fifjuhdn 

JM ^ JM 



< 0, 

where we have used Proposition 4.1 and ht + Ah = 0. 



Claim : 



lim / 1 



lim / vhdfj, < 



Combining (53) and the claim, we have that 

hvdiJ < 



JM 



for any to > and any positive functions h. This implies that v <0. 



□ 



The key point is to prove the claim. Wc need an upper bound on derivatives 
of the logarithm of the heat kernel and the small time asymptotic estimates on 
logarithmic derivatives of the heat kernel. These known results are summarized in 
the Appendix. 



Proof. (Proof of the Claim:) Define 



h{to - t,y)v{t,x,y)dii. 



M 



The proof is along the line of [22]. Wc first show that for any fixed x S M, Wh{t) 
has a finite upper bound as i ^ 0. For the fundamental solution F{t, x, y) = , 
V = [tAf + t{l + kt){Af - |V/n + / - f ln(47ri) - n(l + ^ktf] F. We can write 

Wh{t) = t{2 + kt) I AfFhdn-t{l + kt) I \\/f\^Fhdii+ I (/ - ^ In 47rt) FM/x 

JM JM JM 2 

-n(l + ^ktf I FhdiJL 

2 J M 



= I + II + III + IV 
By Theorem A2 (3) in the Appendix, we have 



VF{t,x,y) 



F{t,x,y) 

AF{t,x,y) 
F{t,x,y) 



dist{x, y) 
t 



+ 



'dist{x, y) 
. t 



+ 
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dist (x,y) , 
n h 



\II\ 



t{\ + kt) J^l^Fhdfi<D\t{l + kt)\Jj 



M 



disl^{x,y) 



Fhdn 



From the asymptotic expansion of the heat hernel of F{t,x,y), and elementary 
computations, we have 



lim 



M 



(54) 



lim / 



M 

< lim 



F{t,x,y)h{y)dii = h{0,x) 
^^ff^F{t,x,y)h{y)df,= ^h{0,x) 



dist{x, y)F{t, x, y)h{y)diJ, 



M 



< lim 

~ t^o 



( /m ^^^^^F{t, X, y)h{y)dt,) (t Fhd/,) = 



where the second limit is by elementary computations on (64), (c.f. pg 8 in [22]). 
Hence, wc have limsup^^Q |/| + \II\ < Dh{0,x). 

When fc > 0, from Cheeger and Yau's heat kernel comparison theorem for Ric > 
—k and Davies and Mandouvalos's lower bound estimate on heat kernel of space 
form, we have 



F{t,x,y) > C(n)(4^i)~"/' 

~ k 

X I 1 + r ' 



expl-;^ 



(n — l)kt (n — X)k 



k 



n — 1 n — 1 
where r = dist{x, y). From F = er^ , wc have 



-t 



l + r 



n — 1 



n, , ^ r"^ (n—l)kt rJ in — l)k , ^, , 
/--ln47rt < - + ^—^ + ^^^—^-lnC{n) 



- ( 



n-1 



Hence, from (54), we have 

lim sup III = lim sup / (/ - - ln(47ri)) Fhdn < - - In C(n) /i(0, x) 

Similar arguments works for the case of fc = 0, where one need to deal with heat ker- 
nel of Euclidean space which has simpler form. On the other hand, limsupj_,o IV = 
—nh{0,x). Hence, we can conclude that limsupf^QWhit) is finite. 

By the entropy monotonicity formula (53), we know that the limit lim^^o W/j(t) = 
7 exists for some finite 7. Hence limt^o ^Vhit) — Wh{^) = 0. By (53) and the 
mean-value theorem, we can find 



such that 
2 



Fhdn = 
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By the Cauchy-Schwartz inequality and the Holder inequality, we have that 
limt.^o ti j (a/ - ^ - ^nk^Fhdfj, 

= ijlimt.^o / i^f - ^^Fhdn- lim ]-nkti [ FhdjjL 

= ti limt^^o / (a/ - ^)Fhdii = 0. 

This yields 

= ti{l-kti) [ {Af-\Vf\^)Fhdfi+ I (f-^\n{4wU) + ^-n{l + htif)Fhdn 



'M JM 

Using integration by parts, we have 



/ (A/ - Vf\')Fhd,i = - f 
Jm J a 

- -j 

Ja 



FAhdfi = -Ah{0,x) 

M 



Hence we have 

lim ti{l + kti) I (a/ - V/ ^)Fhdii = - lim ^^(1 + kti)Ah{0,x) = 
ti— ►o Jj^ \ J ti— >o 

Prom (64) in the Appendix, we have 

hm(/-|ln(4;rt,)-^^^^|^) = - lim^, In ((47rtO"/^e"^F(f„ x, j/)) 

= -lnHo{x,y) 

holds uniformly on any compact subsets of M \ Cut{x), and for y = cxp^(y), 
Hq(x, y) is given by the reciprocal of the square root of the Jacobian of exp^. at Y, 
and Ho{x,x) = 1. Hence we have 



7 = lim / / - - ln(47rti) + - - n(l + -kuf Fhdfi 
ti— >o \ z A z / 

/" / . n, . dist'^{x,y)\ , 

< -{lnHo{x,x))h{0,x) =0 

The last inequality comes from uniformly convergence theorem with Cut{x) zero 
measure for first term, and (54). 

Hence we prove 7 < holds, which is our Claim. 

□ 
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5. Another Harnack Inequality for Heat Kernels 

Combining differential Harnack inequality (19) with the heat equation ut — Au, 
u = n , and |V/|^ — A/ + /j + ^ = we have a Hamilton- Jacobi inequality 

(55) |v/|2 + (2 + A;t)/t + ^ <fc^(2 + H). 

The case of fc = is due to L. Ni [23]. Naturally, the differential Harnack inequality 
(19) will leads to the following Harnack type estimates. 

Theorem 5.1. Let Af" be a Riemannian m,anifold with Ric > —kg (k>0). If we 
denote t := t{t) = 2^^, then for any xi, xi G M", and < < ti, the following 
Harnack type estimates hold, 

(56) Vf2fiX2,t2) - VTlfixuh) < ^.f^l'^y^^l) + - ^(il)), 

where i\ = i{t\), ^2 = ^(^2); 

^{t) = k[ Vtdt = k~i [s/ktikt + 2) - ln(l + kt+ ^/kt{kt + 2))], 
Jo 

and s = s{t) is defined as following 

(1) s{t) := / — = A;-5 [^ykt{kt + 2) +ln(l + kt+^/kt{kt + 2))\, for k > 0; 

Jo vt 

(2) s{t) := y/i, for k = 0. 

In this section, we always assume {M",g) a complete (possibly noncompact) 
manifold with Ricci curvature bounded from below, i.e. Ric{M) > — fc for some 
constant k > 0. We shall apply an equivalent form of our differential Harnack 
inequality (19): 

t{2 + kt)Af - t{l + fci)|V/|^ + / - n(l + ^ktf < 

for the heat kernel F{t, x, y) = „ on M to obtain a pointwisc Harnack inequal- 
ity The estimate (19) is a Li-Yau-Hamilton type since combining with the heat 
equation 

A = A/-|V/|2-^, 
we have a Hamilton-Jacobi inequality 

f n 

(57) \Vff + {2 + kt)ft + ^<k-{2 + kt). 
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Proof. (Proof of Theorem 5.1:) We only prove the case of fc > here. Case 
fc = is due to Ni [23]. Let i{t) and s{t) be defined as in the theorem. 

Vi^fit2,X2,y)-Vi'if{ti,xi,y) = ^ j^(Vif{t,-fit),y)yt 

= fj Vi[ft + ^/ + {vf,^'{t)))dt 

< rsfi(-^\Vff + \Vf\-W\ + '$)dt 



< 



2^1^/1^ +i^/Hyi + ^). 

' (iVi(2 + fci)|7'(i)l' + \nkVTjdt 



for any path -f{t) joining from xi to X2- One can check that s'{t) = This gives 
a Harnack type estimate: 

Vt2f{t2,X2,y)-Viif{tuxi,y) < inf-/ |y(s)|2ds + -(^(ts) - $(ii)) 

Choose 7(s) to be a shortest geodesic with constant speed completes the proof. □ 

Remark 5.1. The k = case was obtained by L. Ni in [23]. The heat kernel 
comparison theorem of Cheeger and Yau [2] will imply the differential Harnack 
inequality as we did in proof of Theorem 1.6. Reversely, when the m,anifold is 
Ricci nonnegative, Ni recovered Cheeger-Yau's heat kernel comparison theorem by 
the above Harnack inequality. We thank Lei Ni for helping us understand his paper 
and pointing out a mistake in the early manuscript. 



Along the line of Ni, we consider the case of A; > 0. We take a;i = o in the above 
theorem, where o is the singular point of the fundamental solution at t = 0. Argue 

Vtf{o,t) 

as in [23], one gets lim vif(t,o,o) — —^(t) < 0, since limu(o,t) = hm , „ = 

i5o(o). This yields the following heat kernel comparison theorem. 

Corollary 5.1. Let M", F, f, t, s, and $ be the same as in Theorem 5.1. Then 
for any {x, y) G M x M, we have 

(1) When Ricci{M) > 0, f{x,t) < ^SSjSI, 

(2) When Ricci{M) > -k, f{x,t) < ^ + f ^. 

6. Entropy formulas with monotonicity 

In this section, we will introduce various new entropy functionals and discuss 
their monotonicity along the heat equation on any compact Riemannian manifold. 
Point-wise differential Harnack inequalities and monotonicity formulas for entropy 
functionals are closely related. Usually, a point-wise differential Harnack quantity 
easily yields a monotonicity formula for the related functional. But reversely, it 
is more difficult. In general, the proofs of monotonicity formulas for functionals 
are also easier. The reason is upon integration over closed manifolds, all the in- 
formation of a divergence form will be disappeared. This implies the point-wise 
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differential Harnack quantity should be the representative in the space of the en- 
tropy integrands for the same entropy functional. 

In this section, based on the Li-Yau type and Li-Yau-Hamilton-Perelman type 
of differential Harnack inequalities we introduced in section 2 and 4, we can easily 
establish monotonicity formulas for the related entropy functionals. But in our 
actual searching for differential Harnack, we discovered the functionals first, then 
localized them and obtained the pointwise version. 

As in section 4, we will follow the notations of Yau [31] and Perelman [24] to 
J. -/ 
assume u = e~' for positive heat solutions and u = — — n- for positive heat kernels. 

^ (47rT)2- 

We introduce the following Li-Yau entropy functional Wly, where the integrand 
is t'^Fu or sinh^{kt)Fu from Proposition 2.1, 
(58) 

yVLY{u, t) := / t^FudjjL = - / ^ In u + — + — (1 + -fcf) udjjL, 



WLY{u,t) := [ smh^{kt)FudiJL = - [ sinh^(t) [Ainu + ^ [coth{kt) + l] 

where we have used integration by parts to get J^^ ftudfi = /^(|Vlnu|^— Alnu)^^/^ : 
0. 

As a direct consequence of Proposition 2.1 and the differential Harnack inequality 
in Theorem 1.1, we have the following theorem. 

Theorem 6.1. Let M" be a closed manifold. Assum,e that u is a positive solution to 
the heat equation (1) with ud^ — 1 and let f ^ — Inn. Consider the functional 



udfi, 



Jm 



Wly = - / t 



» , n nk ,^ 1 , ■ 



udfi. 



IfRicci{M) > -k, then Wly < for allt>0 
(59) 

-rWLY{f,t) = -2e [ \fij-^gij-^^gij\^udti-2t'' [ {Rij+kgij)fifjudii < 0. 

The monotonicity is strict for all t > 0, unless the manifold is Einstein, i.e. 
Ricci{M) = —k, and f satisfies the gradient Ricci soliton equation \Rij + fij — 

2i9ij = 0- 

Remark 6.1. On a complete noncompact manifold, the above Ricci soliton equation 
with potential function as the logarithmic of a positive heat solutionin can be realized 
in some cases, e.g., heat kernels on M" obtains the equality with k = 0. 

Proof Recall that /* = A/ - [V/l^. Using (8), we have 

Wly = / t'^Fud/j. < 0. 
Jm 

The monotonicity follows from Wly = / t^Fudfi and Proposition 2.1. □ 

Jm 
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Exactly the similar theorem is also true for 



Wly {u, t) = - [ sinh^ {kt) [a In u + ^ [ coth(fcf) + l] 



udji. 



We leave the details to the readers. 



On the other hand, one can prove these theorems directly. Since the entropy in- 
tegrand becomes simpler than the differential Harnack quantity and integration by 
parts works for closed manifolds, one can get an easier and more direct derivation 
for the monotonicity formula. We will use this idea to derive the Perelman type of 
entropy monotonicity in the proof of Theorem 6.2. 

In regard of the Perelman type LYH Harnack quantity, the following entropy 
formula is very natural. We define 



(60) Wp(/, r) = j^^ (r| V/p + / - n(l + Ikrf) ^ 

with / -^dn = 1. When fc = 0, Wp is exactly Ni's functional W in [211. The 

Jm (47rT) 2 

following theorem generalizes L. Ni's result in the sense that for closed manifolds 
there is no curvature condition needed. 

Theorem 6.2. Let M" be a closed manifold. Assume that u is a positive solution 
to the heat equation (1) with udfi = 1 . // we choose k E ^ to he any constant 

satisfying Ricci{M) > —k and let f be defined (IS U — -rr and T = T{t) with 

^ = 1, then Wp < for all t > 0, and 
(61) 

Awp = -2t^^^ \i9ij+fi3-l^9i3?j^dtx-2T j^^{Rij+kgij)fifjj^-^ < 0. 

Moreover, the monotonicity is strict for all t > 0, unless the manifold is Einstein, 
Ricci{M) = —k, and f satisfies the Ricci soliton equation ^Rij + fij — ^9ij = 0. 

There are various proofs for entropy monotonicity formula. An immediate proof 

is by using Proposition 4.1. Here, wc will present a direct proof which is based on 
a change of variable argument, see similar argument in [14]. 

Proof (Proof of Theorem 6.2) Let / := -In?! = / + f ln(47rT). Easy to 
see ft — A/ — |V/p. We first observe that the derivative of the Nash entropy 

AA(/,t) := / fe-f'd^l is 
Jm 

^{f,r)= I ftudix+ j futdtJi= j {Af-\Wff)udfx+ [ fAudfx 
Jm Jm Jm Jm 

= / Afud/j,, 
Jm 

where the last step follows from integration by parts. As far as integration by parts 
is allowed, we can use the fact Jj^ \Vf\^e~fdfj, = Jj^ Afe~^diJ,. 



It 
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Applying Bochner formula, one get, 

= / Afudfi + T I (A/(W + Afut)diJ. 
Jm Jm 

= [ Afudn + T [ (2A/- IV/HAud/i 
Jm Jm 

= / AfudiJ,-2T {\fij\^ + Rijfifj)udn 
Jm Jm 

= I Afudn-2T [ \fij\^udii + 2kT [ iV/pwd/x 
Jm Jm Jm 

-2t / {Rij +kgij)fifjudfi 
Jm 

Now we are ready to obtain the monotonicity by using integration by parts and 
completing the square. 



d_ 
It . 



■ (t|V/T + / - I ln(47rr) - n(l + ifcT)2)«d/x 

f f f Tl Tlfc^ 

= 2/ Afudn-2T \fij\'^udii + 2kT \V f\^udix - — + nk — r 

Jm Jm Jm 2t 2 

-2r / {Rij + kgij)fifjUdii 
Jm 

= -2t / \fij - (77- + 7;)9i3?udii -2t I {Rij + kgij)fifjudii 
Jm ^ Jm 

Change / back to /, we complete the proof. 



□ 



The third interesting entropy functional is the 'Nash entropy', — jj^ HlogHd^i, 
where H is the positive heat kernel. Wc will use the linearized version of our 
generalized Li-Yau estimate, namely, the estimate in Theorem 1.2, to illustrate the 
idea. The nonlinear version works exactly in the same way. 

Following the ideas in Section 5 of [24] and motivated by and along the line 
of Addenda to [21], we discuss the relations among these different entropies. Let 
u{x,t) be a positive solution to the heat equation with udfi = 1. We define 



N{u,t) = / —{\ogu)udfi 
Jm 



IM 

and 



~ 71 n \ 71 

N{u,t) = N{u,t) - -log(47rt) - -kt{l + -kt) - -. 
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Direct computations shows that 

-/ (A(logu) + ^ + ^(l + jH))«<iC 
/J|V,og„|^-(l + |«)^-^-f(l + H.*, 

where in the last step we have used integration by parts and the heat equation. 

Notice that the integrand in the last step is just the generalized Li-Yau gradient 

estimate (8), which is 

(63) iviogni^ - (1 + Ikt)^ -Yt-Y^^ + 1'''^ ^ °' 

for any closed manifold when choosing proper k. Using (63), one arrives at the 
following estimate on the Nash entropy — Jj^HXogHdji, which extends L. Ni's 
result to general manifolds. 

Proposition 6.1. Let M" he a complete Riemannian m,anifold with Ricci{M) > 
—k and H be the positive heat kernel. Then N{H,t) satisfies the following proper- 
ties: 

(1) -^N < 0, unless M is an Einstein manifold and H satisfies the gradient 
Ricci soliton equation \Rij — (liiu)y — -^gij = 0. 

(2) limt^oiV(ff,i) =0. 

Proof. The monotonicity is a simple consequence of the generalized Li-Yau gradient 
estimates for heat kernels on complete manifolds with Ricci{M) > —k. Since the 
manifold may be noncompact, one can not use Theorem 1.2 directly. But for heat 
kernels, one can easily extend Theorem 1.2 to noncompact manifold by using the 
techniques we developed to prove the claim in Section 4. The equality case is from 
the vanishing of the first variation formula. 

□ 

The study of relations between pointwisc differential Harnack inequality and 
monotonicity of entropy functionals for Ricci flow equations and heat equations is 
an important and very active field. As we have revealed in this paper, for both 
equations, Ricci soliton plays important role. Sec Entropy formulas for Ricci flow 
in Perelman's original work [24], and others, e.g., [9], [14]. 

7. Appendix 

We summarize some known results about heat kernels on manifolds in this Ap- 
pendix. 

Theorem A.l. Let M" he a smooth, complete Riemannian manifold. Let C{M) c 
M X M he the set of pairs of points {x,y) such that y G Cut{x). Let F{t. x, y) he 
the positive fundamental solution of the heat equation dtu{t,x) = Au{t,x), define 

Et{x,y) = -2tlnF{t,x,y), E{x,y) = ^dist'^{x,y). 



(62) 



dN 
dt 
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Then there are smooth functions Hi{x,y) defined on (M x M) \ C{M) such that 
the asymptotic expansion 

(64) F{t,x,y)^[—\ e ^ ^ i/,(a;, 2/)f 

holds uniformly as t on compact subsets of (M x M) \ C{M). Further, if 
y = exp^{Y), then Ho{x,y) is given by the reciprocal of the square root of the 
Jacobian of exp^ at Y. 

Furthermore, the following estimates on logarithmic derivatives of the heat ker- 
nel on M are known: 



Theorem A.2. 

(1) (Varadhan [29], Cheng-Li- Yau [3]) On any compact subsets of M x M, 

lim Et{x,y) = E{x,y) uniformly; 

(2) (Malliavin and Stroock [19], Stroock and Turetsky [27]) On any 

compact subsets of (M x M)\ C{M), 

limV"'Et{x,y) = V"'E{x,y) uniformly; 

(3) (Stroock and Turetsky [28], Hsu, Elton [13]) There are upper bounds 

for derivatives of the heat kernel on any closed manifold M as 

\V"^Fit,x,y)\ < D^{^^^^^ + l.yFit,x,y) 

where the are some constants depending only on M . 

(4) (Neel [20]) For any A e TyM, we have 

-\A\dist{x,y) < liminft^o V^-Ef (x, y) 

< limsupt^o VA-Et(a;,2/) < \A\dist{x,y) 
and 

\A\'^dist^{x, y) < liminft^o tW\ j^Et{x, y) 

< limsupt^o*V^,A-^t(a;,2/) < 
hold for any {x, y) € M x M. 

For heat kernels on hyperbolic space, Davies and Mandouvalos have the follow- 
ing estimates. 

Theorem A. 3. (Davies-Mandouvalos [8]) Let F^{t,x,y) = F^ {t,dist{x,y)) 
be the heat kernel of A on , the space form with constant sectional curvature 
-K < 0. Then 

c(n)"^/i(t, dist(x,y)) < F^{t,x,y) < c{n)h{t,dist{x,y)) 

where c(n) depends only on dimension n and 

{n-lfKt {n-l)^/Kr\ 



(65) h{t,r) = (47rt)-"/2exp(-^- 



x(l + y/Kr + Kt^ ' (l + y/Kr). 
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Another very useful estimate is the heat kernel comparison theorem of Cheeger 
and Yau. 



Theorem A.4. (Cheeger- Yau [2]) 

F{t, X, y) > F^{t, x,y),y Ric{M) >-{n- 1)K. 
We have the following lower bound estimates for heat kernel on M: 

(66) F{t, X, y) > c{n)~'^h{t, dvst{x, y)) 

where h{t,r) as in (65). 

The following theorem is from [26] page 167, Corollary 2. 

Theorem A. 5. (Li- Yau [26]) Let H{x,y,t) be the heat kernel of a complete 

Riemannian manifold M. For any p > 0, T > 0, set 

(67) F{y,t)= [ H{x,^,T)H{^,y,t)d^. 

Jm\b^(p) 

Then for any 5 > Q, and R> 0, 



(68) -^^-(^^ 



/ F^{y,il + S)T)dy 
< exp (^) ■ exp (2(1^) • ^ ^ H\x,i,T)d^. 



lM\B^{p) 

Moreover, if p = 0, i.e. 

(69) F{y,t)= [ H{x,^,T)H{^,y,t)d^, 

JM 

then for any (5 > 0, T > 0, and R > 0, we have 

(70) ^ F^y, (1 + 6)T)dy < exp (^)i?(x,T). 
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